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1.  INTRODUCTION 

Consider  the  standard  multivariate  linear  regression  model  (cf.  Anderson 
(1984),  Chapter  8): 

Z  =  (5£>  +e,  (1.1) 

where  Z  :  p  xN  is  the  matrix  of  observations,  (3 :  p  xq  is  the  matrix  of  unknown 
regression  coefficients,  D  :  q  xN  is  the  design  matrix,  and  e:  p  xN  is  the  matrix  of 
unobservable  random  errors.  Assume  that  q  <N,  D  is  of  full  rank  q ,  and  that 


e  -  N(0,  Z®/*),  (1.2) 
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which  indicates  that  the  N  columns  of  e  are  mutually  independent  p  -variate  normal 
random  vectors  with  zero  mean  and  common  unknown  covariance  matrix  I,:  pxp, 
assumed  positive  definite.  We  shall  study  the  consistency  of  invariant  tests  of  the 
general  linear  hypothesis 

H0:  pD0  =  0,  (1.3) 

where  D 0:  q  xr  has  full  rank  r  (1  <r  <q).  If  p  is  partitioned  as  (P1(  P2)  with 
Px :  q  xr,  an  important  special  case  of  (1.3)  is  H 0:  px=0. 

The  general  multivariate  analysis  of  variance  (MANOVA)  testing  problem  is 
that  of  testing  (1.3)  vs.  (1.1).  It  is  well  known  that  this  testing  problem  can  be 
reduced  by  sufficiency  and  invariance  to  the  following  canonical  form  (cf.  Anderson 
(1984),  Section  8.3.3  or  Lehmann  (1986),  Sections  8.1,  8.2):  based  on  the 
independent  observations 

X(pxr)  —  N(£,L  ®/r) 

Y(pxn)  —  N(0,Z®/„),  (1‘4) 

where  § :  p  x  r  is  a  matrix  of  unknown  means,  test 

H0:  £=0  vs.  H:  (Z  unknown).  (1.5) 

-  1 

We  assume  that  n  (=N-q)>p,  so  that  L=  —  YY  is  positive  definite  with 

n 

probability  one. 

The  canonical  testing  problem  given  by  (1.4)  and  (1.5)  is  invariant  under  the 
group  of  nonsingular  linear  transformations 

(X,Y)  (BXV^BY' ¥2),  (1.6) 

where  B  :  pxp  is  nonsingular  and  r  xr,  XP2:  nxn  are  orthogonal.  Under 
(1.6),  the  parameters  of  the  model  (1.4)  are  transformed  according  to 

&£)  (B^.BZB').  (1.7) 

The  maximal  invariant  statistic  and  maximal  invariant  parameter  may  be 
represented  as 
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We  shall  be  concerned  with  the  consistency  of  invariant  tests  for  (1.5)  =  (1.13),  i.e., 
tests  that  depend  upon  ( X,Y )  only  through  c  (or,  equivalently,  through  d)  and 

whose  power  functions  therefore  depend  upon  (£,  X)  only  through  X  =  (Xj . X,), 

the  vector  of  noncentrality  parameters.  Since  c,  estimates  X,,  a  “good”  invariant 
test  should  accept  (reject)  H0  for  small  (large)  values  of  c  j. . . . ,  c,  (equivalently,  of 
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dv  ...  ,dt).  In  fact,  Schwartz  (1967b)  has  shown  that  every  admissible  invariant 
test  for  (1.5)  =  (1.13)  must  have  a  monotone  acceptance  region  A  in  terms  of  c  or 
(equivalently)  d.  That  is  (in  terms  of  d),  if  d  =  (dlt  ...,dt)  e  A  cD,  and 
d,m(dl',...,dt")e  D,  is  such  that  d'<d  (i.e.,  dx'<  dx, . . .  ,d,'<dt),  then 
d'  e  A.  Therefore,  we  shall  restrict  our  attention  to  the  class  of  monotone  invariant 
tests,  i.e.,  those  with  monotone  acceptance  regions. 

Perlman  and  Olkin  (1980)  showed  that  every  monotone  invariant  test  is 
unbiased  for  testing  H0  vs.  H.  The  criterion  of  unbiasedness,  therefore,  does  not 
distinguish  among  admissible  invariant  tests.  Likewise,  neither  does  the  classical 
notion  of  consistency,  which  we  shall  call  sample  size  consistency  (SSC).  In  this 
paper  we  introduce  the  notion  of  parameter  consistency  and  show  that  it  does 
distinguish  among  admissible  invariant  tests. 

An  invariant  level  a  test  with  acceptance  region  A  is  said  to  be  parameter 
consistent  (PC)  if,  for  fixed p,r,n,  and  a,  its  power 

Px,{ddA}  -»  1 

as  one  or  more  noncentrality  parameters  Xt  — >  <*>.  It  will  be  seen  that  the  well- 
known  Bartlett-Nanda-Pillai  trace  test,  which  is  both  admissible  and  the  locally  most 
powerful  invariant  test  for  H0  vs.  H,  fails  to  be  PC  unless  a  or  n  is  sufficiently 
large,  whereas  the  Roy  maximum  root  test,  the  Lawley-Hotelling  trace  test,  and  the 
likelihood  ratio  test  (=  Wilks  criterion)  are  both  admissible  and  PC  for  every  a  and 
n  (cf.  Section  4). 

It  is  important  to  notice  that  parameter  consistency  is  defined  in  terms  of  the 
power  of  a  single  invariant  level  a  acceptance  region  A  at  sequences  of  alternatives 
{X}  with  ||7.||  — >  oo,  whereas  sample  size  consistency  is  defined  in  terms  of  the 
limiting  power  of  a  sequence  of  invariant  level  a  acceptance  regions  {A(n)}  at  a 

fixed  alternative  X*  *  (0 . 0).  Usually  the  sequence  {A is  defined  in  terms  of 

a  single  invariant  test  statistic  /  =f  (d)  as  follows: 

A/">  =  Af(Ca)  *  {deDt\f(d)<ca},  (1.14) 
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where  ca  =  ca(p,r,n;f)  satisfies 

Px^o{f(d)<ca}  =  1-a.  (1.15) 

If  /  is  monotone  on  D,  (i.e.,  nondecreasing  in  each  dt)  then  A_r(ca)  is  a  monotone 
invariant  acceptance  region  with  power  function  given  by  P\{f  (d)  >  ca}. 

Necessary  and  sufficient  conditions  for  the  parameter  consistency  of 
monotone  invariant  tests  are  presented  in  Section  2,  while  sample  size  consistency  is 
defined  and  characterized  in  Section  3.  The  relation  between  parameter  consistency 
and  sample  size  consistency  of  monotone  invariant  tests,  in  particular  admissible 
invariant  tests,  is  examined  in  Section  4.  Few  detailed  proofs  are  given,  as  this 
paper  is  primarily  expository.  The  proofs,  together  with  extensions  of  the  results  to 
related  multivariate  testing  problems,  will  appear  in  Anderson  and  Perlman  (1988). 

2.  PARAMETER  CONSISTENCY  OF  MONOTONE  INVARIANT  TESTS 

In  a  general  hypothesis-testing  problem,  a  level  a  test  of  H0  vs.  H  is  said  to 
be  parameter  consistent  if,  for  fixed  sample  size,  its  power  approaches  one  for 
sequences  of  alternatives  in  H  whose  Kullback-Leibler  discrimination  distance  from 
H0  becomes  arbitrarily  large.  For  the  canonical  MANOVA  problem  (1.5)  =  (1.13), 
this  definition  is  equivalent  to  the  following: 

DEFINITION  2.1.  For  fixed  p,  r,  n,  and  a,  an  invariant  level  a  test  for  (1.5)  ■ 
(1.13)  is  parameter  consistent  (PC)  if  its  power  at  A  =  (Aj, . . .  ,  A,)  approaches  1  as 
||A||  where  ||A||  =  Z/A;  =  tr^'Z-1.  For  i  =  1, ...  ,t,  the  test  is  parameter 
consistent  of  degree  i  (PC  (i ))  if  its  power  at  A  approaches  1  as  A,-  —>  °°.  □ 

Since  X>  £A,5  0,  obviously  PC(/)=> PC(i+l),  and  PC<s=s>PC(l). 
It  will  be  seen  in  Section  4  that  parameter  consistency  is  not  equivalent  to  sample 
size  consistency. 

In  order  to  study  the  power  of  an  invariant  test  at  the  alternative 
X=  (Xj . A,)  we  may  assume  that  Z)  =  (n,Ip),  where  p: p  xr  is  any  matrix 
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such  that  chi(\ 41')  =  X;,  1  <i  </,  and  where  Ip  is  the  p  xp  identity  matrix.  Under 
this  assumption  XX '  and  YY'  have  standard  Wishart  distributions,  noncentral  and 
central  respectively,  and  are  independent.  Define 

/,  =  lt(X)  =  chAXX'),  1  <i<t. 

Our  characterization  in  Theorem  2.3  of  parameter  consistency  for  monotone 
invariant  tests  is  based  upon  the  following  technical  result: 

LEMMA  2.2. 

p 

(i)  /,•  00  iff  X,-  — » 

p 

(ii)  c(-  — »  00  iff  X,-  — »  00. 

p 

(iii)  di  1  iff  X,-  — >  00. 

PROOF.  The  result  (i)  follows  from  appropriate  stochastic  bounds  for  l-t  in  terms  of 
Xj,  (ii)  follows  from  (i)  by  conditioning  on  Y,  while  (ii)  and  (iii)  are  equivalent  by 
(1.1 1).  See  Anderson  and  Perlman  (1988)  for  details.  □ 

It  is  most  convenient  to  state  our  results  for  acceptance  regions  A  defined  in 
terms  of  the  statistic  d.  For  any  subset  A  c  D,  (recall  (1.12)),  we  denote  the  closure 
of  A  in  D,  by  A ,  where 

D,  =  {x  |  l>xj>  •••  £x(>0}.  (2.1) 

If  A  is  monotone  then  A  is  also  monotone  and  A  \  A  has  Lebesgue  measure  0. 
Since  the  distribution  of  d  is  absolutely  continuous  with  respect  to  Lebesgue 
measure  for  every  Xe  H0\jH,  the  power  (and  hence  the  consistency)  of  the 
invariant  test  with  acceptance  region  A  is  the  same  as  that  of  the  test  with 
acceptance  region  A . 

Let  e  0,  e  . . e,  denote  the  vertices  of  D„  i.e.. 


i  t-i 


and,  for  0<T|  <  1,  define  e,(T|)  e  D,  by 
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(2.3) 


<?,(*!)  =  «r +*!(*, -.«,•) 


i  l-i 


If  A  is  a  monotone  subset  of  D„  then  e,(T|)  e  A  implies  «,  _i(T|)  e  A . 


THEOREM  2.3.  Fix  p,r,n,  and  a,  let  A  cD,  be  a  monotone  level  a  acceptance 
region  for  the  testing  problem  (1.5)  a  (1.13),  and  fix  i  e  (1, . . . ,/}.  A  necessary 
and  sufficient  condition  that  the  invariant  test  with  acceptance  region  A  be  PC  ( i )  is 
that  (Ti)  4  A  for  all  T|  >  0.  A  sufficient  condition  is  that  ei<£  A.  □ 

If  the  upper  boundary  of  A  is  not  too  irregular.  Theorem  2.3  essentially  states 
that  the  test  with  acceptance  region  A  is  PC(i)  if  either  eL  4  A  (cf.  Fig.  2.1a)  or  e, 
lies  in  the  upper  boundary  of  A  (Fig.  2.1b),  whereas  it  fails  to  be  PC  (i)  if  e,  e  A 
but  4  (upper  boundary  of  A)  (Fig.  2.1c).  Thus,  the  test  is  PC (i )  if  e;  lies  either 
in  the  rejection  region  or  in  its  lower  boundary.  These  three  cases  are  illustrated  in 
Figures  2.1a,b,c  where  i  =  1  and  t  =2. 


Fig.  2.1a: 
testisPC(l) 


Fig.  2.1b:  Fig.  2.1c: 

testisPC(l)  testis /wtPC(l) 


Note  that  if  a  >  0,  every  monotone  invariant  test  must  be  PC  (r ). 

It  is  convenient  to  restate  Theorem  2.3  for  the  case  where  A  =Af{ca)  as  in 
(1.14)  and  /  is  a  monotone  invariant  test  statistic  defined  on  Dt.  For  any  extended 
real-valued  function  /  on  D,  define  /  on  D,  as  follows: 

/(*)  =  f(x~)  =  lim/((l-e)x) .  (2.4) 

e-lo 

Then  /  is  defined  (possibly  extended  real-valued),  monotone,  and  lower 
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semicontinuous  on  D„  hence  Aj{c„)  =  Af{c^)  is  a  closed  and  monotone  subset  of 
D,.  Furthermore,  {x  a  D,  \  f  (x)*f  (x)}  has  Lebesgue  measure  zero,  hence  so 
does  Af(c^)\  Ay(ca)  for  every  cw  Thus,  the  level  a  tests  determined  by  /  and / 
are  equivalent  and  have  identical  power  functions. 

COROLLARY  2.4.  Fix p,r,n,  and  ct,  let /  be  a  monotone  test  statistic  defined  on 
D„  let  ca=  ca(p,r,n;f)  satisfy  (1.15),  and  fix  /  e 

(i)  A  necessary  and  sufficient  condition  that  the  invariant  level  a  test  based  on  /  be 
PC  (/ )  is  that  /  (c;  (T|))  >  c  a  for  all  T|  >  0.  A  sufficient  condition  is  that  /  («,■ )  >  c  a. 

(ii)  If  1)  =  0  is  a  point  of  increase  of  /(e;(T|))  (i.e.,  /  (e,  (q))  >  /  (e,)  for  all  tj  >  0) 
then  /  (e,)  >  c  a  is  a  sufficient  condition  that  the  level  a  test  based  on  /  be  PC  (/ ).  If 
/  («»(n))  is  continuous  at  Ti=0  then  /  (e,)  >  ca  is  a  necessary  condition.  Thus,  if 
T) =0  is  both  a  point  of  increase  and  a  continuity  point  of  /  (e;(T))),  then  /(e;)  >ca  is 
both  necessary  and  sufficient  for  the  level  a  test  to  be  PC  (/ ).  □ 

If  /(d)  depends  only  on  di+1, . . .  ,dk  for  some  0<i  <k  <t,  it  follows  from 
Corollary  2.4(i)  that  for  0<a<l,  the  level  a  test  based  on  /  is  PC(£)  but  not 
PC  (/).  It  may  or  may  not  be  PC  (J)  for  i  <j  <k:  for  example,  the  level  a  test 
based  on  f(d)  =  *s  PC  (A)  but  not  PC(A-l),  while  those  based  on 

UL^-djT1  and  2*+i dj(\-d.j)~x  are  PC (i  + 1)  but  not  PC (i ).  Furthermore, 
parameter  consistency  may  depend  on  the  value  of  a:  the  level  a  test  based  on 
/(^)  =  Xi+t^y  k  PC  O')  but  not  PCO'-D  for  a  satisfying  j-i-1  <ca<  j-i, 
j  =  i+ 1, . . . ,  k .  Further  examples  are  considered  in  Section  4. 

3.  SAMPLE  SIZE  CONSISTENCY  OF  MONOTONE  INVARIANT  TESTS 

For  fixed  p,  r,  and  a,  we  shall  study  the  consistency  of  the  sequence  of 
invariant  level  a  tests  determined  by  the  acceptance  regions  Ajn\n  >p)  based  on  a 
monotone  test  statistic/  =  /  (d)  defined  on  D,  (see  (1.14)). 


DEFINITION  3.1.  For  the  testing  problem  (1.5)  =  (1.13),  the  sequence  of  invariant 
level  a  tests  based  on  /  is  sample  size  consistent  (SSC)  at  the  fixed  alternative 
X*  *  (0, ....  0)  if  the  power 

P \M{f  (d)  >  ca(p,r ,n;f)}  1 

as  n  — >  oo  for  every  sequence  of  alternatives  {X(n)}  such  that  X(b)  =  (n+o(n  ))X*, 
i.e.,  n-1X(n)  — >  X*.  If  this  holds  for  every  0<a<l,  we  say  that  /  is  sample  size 
consistent  at  X*  (SSC  at  X*).  The  test  statistic  f  is  called  sample  size  consistent 
(SSC)  if  it  is  SSC  at  every  X*  *  (0 . 0).  □ 

REMARK  3.2.  The  condition  X(n)  =  (n  +o(n  ))X*  stems  from  the  fact  that  X  in 
(1.4)  is  typically  of  the  form  ■'//T(l+o(l))X  with  X  a  sample  mean  vector.  In  the 
simplest  case,  for  example,  X\, . . .  ,XN  are  independent  univariate  observations 
from  N((x,  a2)  and  we  wish  to  test  |i=0  vs.  |i*0.  Here  the  two-sided  t  -test  rejects 
|X= 0  for  large  values  of  T2  =  (Jn  X)2/ s2,  where  s2 ~'L{Xi-X)2 ! N In  this 
case,  p=q-r-\,  t=l,  n  =N- 1,  and,  when  p*0,  T2  has  the  noncentral  F  distribution 
with  1  and  n  degrees  of  freedom  and  noncentrality  parameter  X  =  N  i^/o2.  This  is  of 
the  form  X=  (n +o(«))X*  with  X*  =  □ 

Theorem  3.5,  our  main  result  on  the  sample  size  consistency  of  a  monotone 
invariant  test  statistic  / ,  is  based  on  the  following  two  elementary  lemmas,  which 
summarize  the  limiting  behavior  of  d  and  ca  =  ca(p,r,n;f)  as  n  with  p,  r, 
and  a  fixed.  Their  proofs  follow  directly  from  the  definitions  (1.4)  and  (1.11)  of  dit 
X,  and  Y. 

LEMMA 3.3.  Fixp,r,X*e  C„andi  6  (1 . f}. 

(i)  If  X= 0,  then  dt  -  Op(l/n )  as  n  — » 

p 

(ii)  If  X  =  («  +o(n ))X*,  then  d  — »  8*  as  n  -»«>,  where 

5*s8(X*)s(Sr . 8,*)  e  D, , 

5f=xf/(xr+D  □.  (3'1) 
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ForO<Ti<l,  define 

/oCn )  ■  /(*o(i 1))  ■  /(n . 11)  •  (3.2) 

To  avoid  technicalities,  we  shall  assume  that 

/0  is  continuous  and  stricdy  increasing,  (3.3) 

hence  the  inverse  function  /0-1  is  well-defined,  continuous,  and  stricdy  increasing. 
Since 


fo  (4)  ^  fid)  <  f0(di)  (3.4) 

ford  s  (rfj, ...  ,d,)  e  Dr,  it  follows  from  (1.15)  that 

^x=oW>/o-1(ca)}  <  a  <  /’x=o{^i>/o"1(ca)}  .  (3.5) 

Because  Fx=0{0<df  <d1cl}=  1,  this  implies  that  0  < /o-1(Ca)  <1  for  0<a<  1, 
which,  together  with  Lemma  3.3(i),  yields  the  following  result; 


LEMMA 3.4.  Fixp,r  anda(0<a<l).  Then 

0  <  liminfn  /0-1  (c a)  <  limsup  n  /0_1  (c „)  <  »  .  □ 


For  8  s  (Sj, ...  ,8, )  €  Dt,  define 


/(§)  =  Urn  /o"1  [/  (8  ve  0(rj))]  (3.6) 

■nio  t) 

provided  the  limit  exists  (possibly  infinite),  where  x  Vy  =  (xjVyj, . . .  ,x,vyt). 
Then  by  the  monotonicity  of  / , 

/  is  monotone  on  D, ,  (3.7) 

1  =  /( 0)  <  /(8)  <  co  (3.8) 


By  (3.7), 


f\b~)  <  fib)  <fib+), 


where  /(8-)  is  defined  as  in  (2.4)  and  where 


(3.9) 
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/(8+)  =  lim/((l+£)8). 
eio 

By  (3.7),  f  (5±)  exists  provided  that  / ((l±e)8)  exists  for  sufficiently  small  e>0. 
We  say  that/  is  radially  continuous  at  8  if 

/(S-)  =/(S)  =/(8+).  (3.10) 

Radial  continuity  is  a  weaker  requirement  than  continuity:  for  example,  if  /’ 
depends  on  8  =  (81( . . . ,  8f)  only  through 

rank  (8)  =  number  of  nonzero  8, 

=  max  {/  |  8,  >0}  ,  (3-11) 

then/"  is  radially  continuous  on  Dt  but  not  necessarily  continuous. 

The  following  characterization  of  the  sample  size  consistency  of  /  in  terms  of 
/  may  be  proved  by  applying  Lemma  3.3(ii)  and  Lemma  3.4.  It  also  follows  from  a 
slightly  stronger  result  in  Anderson  and  Perlman  (1988). 

THEOREM  3.5.  Fix  p,  r,  and  X*€  C„  and  set  8*  =  8(X*)  as  in  (3.1).  Let  /  be  a 
monotone  invariant  test  statistic  defined  on  Dt  and  satisfying  (3.3). 

(i)  If  / (8*-)  =  oo,  then  /  is  SSC  at  X*. 

(ii)  If/ (8*+)  <  oo,  then/  is  not  SSC  at  X*. 

(iii)  Suppose,  in  addition,  that  /  exists  and  is  radially  continuous  at  8*.  Then  /  is 
SSC  at  X*  iff  / (8*)  =  oo.  □ 

It  is  important  to  note  that  when  Theorem  3.5  is  applicable,  the  sample  size 
consistency  of  the  sequence  of  level  a  tests  based  on  /  does  not  depend  on  the  value 
of  a  (0<a<  1). 

REMARK  3.6.  In  addition  to  the  hypotheses  of  Theorem  3.5,  suppose  that/0(0)  =  0 
and  that 

/oOl)  ~  <*T\b  as  T|  4. 0  (3.12) 

for  some  a,b  >0.  Then  it  may  be  shown  that  for  any  8  e  D„ 
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(3.13) 


/(S)  =  lim 

TJ-lO 


/(SVe0(T1)) 

foW 


If  (8)] 4  • 


That  is,  / (8)  exists  iff  / (5)  exists,  in  which  case  (3.13)  holds.  Then  all  conclusions 

~  a 

of  Theorem  3.5  remain  valid  with  /  replaced  by  / ,  which  is  usually  easier  to 

calculate.  The  condition  (3.12)  is  satisfied,  for  example,  if  /(d, . dt)  is 

monotone  on  Dt,  and  admits  a  power  series  expansion  about  (0, . . . ,  0).  □ 


LEMMA  3.7.  Suppose  that  /  is  monotone  on  D,  and  satisfies  (3.3).  Then  for 
5e  D„ 


/  (S)  >/o(0)  =>  /(8)  =  ~. 


(3.14) 


PROOF.  It  follows  from  (3.3)  and  the  monotonicity  of  /  that  for  every  M  >0  and 
sufficiendy  small  i)  >  0, 

/(8)>/(0)  ==>  / (8)  >f0(T[M) 

=*>  ~faX  [/(5Ve0(T|))]  >  ^/(f1  If  (5)]  >M  .  (3.15) 

T)  T1 

Now  let  t|  4-  0  and  M  t  oo.  □ 


If  X*  =  X(4*.  L*)  and  5*  =  8(X*),  then  from  (1.8)-(1.11)  and  (3.1), 

rank  (X*)  =  rank  (8*)  =  rank  (£*) .  (3.16) 

COROLLARY  3.8.  Suppose  that  /  is  monotone  on  Dt  and  satisfies  (3.3).  Fix 
X*  e  C,  and  set  8*  =  5(X*). 

(i)  If  /  (8*-)  >  /0  (0),  then  /  is  SSC  at  X*. 

(ii)  If  rank  (8*)  =  t ,  then  /  is  SSC  at  X*. 

PROOF,  (i)  For  sufficiendy  small  e  >  0, 
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/  (5*-)  >/0(0)  =>  /  ((1-  e)5*)  >  /0  (0) 

=^/((l-e)8*)  =  ~ 

by  Lemma  3.7,  hence  /  (8*-)  =  By  Theorem  3.5(i),  /  is  SSC  at  X*. 

(ii)  For  sufficiently  small  e  >0, 

rank  (8*)  =  t  <=o>  rank  ((1  -e)5*)  =  t 
<=£>  (l-e)8(*  >  0 

=>  /  (( 1  -  e)8*)  >  /0((1  -  e)8*)  > /0(0) 

by  (3.4)  and  (3.3).  As  in  the  proof  of  part  (i)  it  follows  that  / (8*-)  =  hence  that 

/is  SSC  at  X*.  □ 

REMARK  3.9.  The  converse  of  Corollary  3.8(i)  is  not  necessarily  true;  in  fact,  it  is 
not  necessarily  true  that  if  /  (8*)  =  /o(0)  then  /  fails  to  be  SSC  at  X*.  For  example, 
if  /  (d )  =  ]"J{  dt  then  /  (8*) = f0  (0)  =  0  whenever  rank  (8*)  <  t ;  however,  f  (8*)  = 

/  (8*)  =  °°  for  every  8*  *  (0 . 0),  hence  /  is  SSC  at  every  X*  *  (0, ....  0)  by 

Theorem  3 .5(i)  or  Remark  3.6.  □ 

The  following  definition  is  suggested  by  Corollary  3.8(ii): 

DEFINITION  3.10.  For  i  =  1... .  ,t,  the  invariant  test  statistic  /  is  said  to  be 
sample  size  consistent  of  degree  i  (SSC(i))  for  the  testing  problem  (1.5)  =  (1.13)  if 
it  is  SSC  at  every  X*  such  that  rank  (X*)  >  i.  □ 

Clearly,  SSC(/)=>  SSC(i+l)  and  SSC  <=>  SSC(l).  If  /  is  monotone  on 
D(,  satisfies  (3.3),  and  is  SSC  at  every  X*  such  that  rank  (8*)  =  i,  then  by  Theorem 
3.5(ii),  / (8*+)  =  oo  for  every  8*  e  D,  of  rank  i.  By  the  monotonicity  of  /  this 
implies  that  /  (8*-)  =  «  for  every  8*  of  rank  i ,  hence  that  f  (8*-)  =  °°  for  every  8* 
such  that  rank  (8*)  >  i ,  and  therefore,  by  Theorem  3.5(i),  that  /  is  SSC  (i ). 
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The  following  result  is  similar  to  Corollary  3.8.  Note,  however,  that  in  part  (i) 
only  the  weaker  condition  /  (8)  >  /0(0)  need  be  assumed,  rather  than  /( 8-)  >  /0(0). 

COROLLARY  3.1 1.  Suppose  that/  is  monotone  on  D,  and  satisfies  (3.3). 

(i)  If  /  (8)  >  /0  (0)  for  every  8  such  that  rank  (8)  =  i ,  then  /  is  SSC  (i ). 

(ii)  If f  {d)  =  g(dx, . . .  ,</,)  for  some  g ,  then /  is  SSC(i)- 

(iii) If  f(d)  =  g(dM,...,dt  )  for  some  g,  then  /  is  not  SSC(i).  In  fact,  /  is  not 
SSC  at  any  X*  such  that  rank  (X9)  <  i. 

(iv)  /  is  SSC(r). 

PROOF,  (i)  By  Lemma  3.7,  /  (8)  =  °°  for  every  8  of  rank  i,  hence  /  (8-)  =  «=  for 
every  8  such  that  rank  (8)  >  i ,  so  the  result  follows  from  Theorem  3.5(i). 

(ii)  If  rank (5)  =  i ,  then  8,-  >  0  and  /(8)  =  gfSj, . . .  ,8,)  >  g(8, . 8.)  = 

/0(8,)  >  /0(0)  by  (3.3),  so  the  result  follows  from  part  (i). 

(iii)  If  rank (8)  <  i,  then  8,+1  =  •  •  ■  =  8,  =  0  and /(Sv*0Cq))  =  g(n, . . .  .tj)  = 
/o(Tl)»  hence  /(8)=1  by  (3.6).  Thus  /(8*+)  =  l  whenever  rank(8*)<i,  so  / 
cannot  be  SSC  at  the  corresponding  X*. 

(iv)  This  is  immediate  from  Corollary  3.8(ii).  □ 

By  Remark  3.9,  however,  it  is  quite  possible  that  /  is  SSC(i)  even  though 
/(8)  =/o(0)  for  every  8  of  rank  i.  To  illustrate  this  more  fully,  consider  the  four 
test  statistics  /  ( d )  appearing  in  the  final  paragraph  of  Section  2.  Each  is  monotone 
on  D„  satisfies  (3.3),  and,  by  Corollary  3.11(ii)  and  (iii),  is  SSC(k)  but  not  SSC(i). 
In  fact,  however.  Theorem  3.5  or  Remark  3.6  implies  that  each  is  SSC  (i+1)  but  not 
SSC(i),  even  though  the  first  two  statistics  311(1  satisfy 

/  (8)  =  /o(0)  for  every  8  of  rank  <  k . 
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4.  COMPARISON  OF  PARAMETER  CONSISTENCY 
AND  SAMPLE  SIZE  CONSISTENCY 

When  comparing  these  two  properties  for  a  monotone  invariant  test  statistic  / 
defined  on  D,,  it  is  important  first  to  examine  their  differences.  Throughout  this 
discussion  the  dimensions  p  and  r  remain  fixed,  while  we  write  d=d(n), 
ca  =  ca(n ),  and  Ay  (c  a)  =  Ay(c  a(n ))  to  stress  the  dependence  of  these  quantities  on 
n  (the  number  of  degrees  of  freedom  for  estimating  £)  —  cf.  (1.4),  (1.11),  (1.14), 
(1.15). 

First  (recall  Definition  2.1),  parameter  consistency  is  defined  for  the  single 
level  a  acceptance  region  Ay  (ca(n))  with  n  and  a  fixed:  we  say  that  Ay  (c„(n))  is 
PC  (i )  (or  simply  “/  is  PC  (i )  for  n ,  a”)  if 

lim  Px{d{n)4  Af(ca(n ))}  =  1 .  (4.1) 

A,— 

It  was  seen  at  the  end  of  Section  2  that  the  PC  ( i )  property  may  depend  non-trivially 
on  the  value  of  cl  By  Corollary  2.4,  this  property  is  determined  not  only  by  the 
values  of  f  {ei  (t]))  but  also  by  that  of  c  a(n )  *  c  a(p ,  r ,  n ;/ ),  therefore  by  the  global 
behavior  of/  onD,. 

On  the  other  hand  (recall  Definitions  3.1  and  3.10),  sample  size  consistency  is 
defined  for  the  sequence  of  acceptance  regions  (Ay(c„(n))  |  n  >p}  for  a  fixed  a: 
this  sequence  is  said  to  be  SSC  (/)  if 

fim  Px»{d(n)4Af(ca(n))}  =  1  (4.2) 

for  every  sequence  of  the  form  A.(b)  =  (n +o(n))X.*  with  rank  (X.*)  >  i . 

Whenever  Theorem  3.5  applies,  this  property  does  not  depend  on  the  value  of  a 
(0<a<l),  so  we  then  simply  say  that  /  is  SSC(/).  Again  by  Theorem  3.5,  this 
property  is  determined  by  the  values  of  / (5)  for  every  5  e  D,  of  rank  i ,  hence  only 
by  the  local  behavior  of  /  in  a  neighborhood  of  the  set  {x  e  D,  |  x,=0}.  (This  is 
because  /  is  always  SSC (t )  (cf.  Corollary  3.11(iv)),  while  for  1<i</-1,  if 
rank  (8)  =  i  then  the  value  of  /  (5)  is  determined  by  the  values  of  /  in  a 
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neighborhood  of  {x  e  D,  |  x,  =  0}. 

In  view  of  these  differences,  it  is  not  surprising  that  the  properties  PC  (i )  and 
SSC(i)  are  not  equivalent,  even  for  a  monotone  test  statistic  /.  In  general,  neither 
property  implies  the  other,  as  demonstrated  by  the  following  examples.  Define  (cf. 
(1.11)) 

fi(d)  =  d,  =  + 

f2(d)  =  'Zdjd-dj)-1  s  trtXY'OTO-1] 
y=i 

h{d)  =  na-rfy)'1  =  detcxx'+yyo/detcyyo 

j= 1 

hid)  =  hdj(l-dj)-1  =  det  (XX  0  /  det  (IT  0 
/=i 

hid)  =  idj  =  tr[XY/(XX’'+yy,)_1] 

y=i 

feid)  =  n  dj  *  det (XX ') /  det (XX '+YY') 

/=i 

Md)  =  d.nd-d/1 
/= 1 

hid)  =  d,  fi(l+dy) 

/= 1 

/9(d)  =  d,  s  cA^JXX^XX'+ITr1]  . 

The  statistics  f\,  h>h>  2nd  /j  are  well-known  (cf.  Anderson  (1984),  Chapter  8): 
/j  is  the  Roy  maximum  root  statistic, /2  is  the  Lawley-Hotelling  trace  statistic, /3  is 
the  likelihood  ratio  (LR)  statistic  (=  Wilks  statistic),  and  f5  is  the  Bartlett-Nanda- 
Pillai  trace  statistic.  Each  of  the  statistics  f\—fg  is  monotone  on  D,  and  satisfies 
(3.3),  (3.10),  and  (3.12)  (more  precisely,  /3  - 1  satisfies  (3.12)).  The  parameter 
consistency  and  sample  size  consistency  of  the  invariant  level  a  tests  based  on 
/i  ~h  are  readily  determined  from  the  results  of  Sections  2  and  3: 
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f  i>  f 2*  fi >  f 4  are  PC(1)  for  all  n ,  a 

and 

f \ >  f 2*  f 2*  f 4  areSSC(l) 

fs  isPC(i),  notPC(i-l),  for 
i— 1  <  ca(n)<  r,i  =  1,  —  ,t 

but 

fs  isSSC(l) 

f6  is  PC(t ),  not  PC(r-l),  for  all  n ,  a 

but 

f6  is  SSC(l) 

/7  is  PC(1)  for  all  n ,  a 

but 

/7  is  SSC(t),  not  SSC(/-1) 

/g,/9  arePC(0. notPC(r-l), forall/i,a 

and 

/g,/9  are  SSC(t ), 

not  SSC(f-l). 


Although  “/  is  PC(i)  for  all  n,a”  and  “/  is  SSC(i')”  are  thus  seen  to  be 
inequivalent,  their  defining  properties  (4.1)  and  (4.2)  have  an  important  common 
feature:  since  X/n)  =  (n  +o(/t))X*,  we  see  that  rank  (X*)  >  i  iff  X/"*  ->  °°  as 
n  ->  oo.  This  suggests  the  following  definition: 

DEFINITION  4.1.  For  fixed  a,  the  sequence  {Af(c a(n ))  |  n>p}  of  level  a 
acceptance  regions  based  on  /  is  said  to  be  eventually  parameter  consistent  of 
degree  i  ( eventually  PC(i))  if  there  exists  n0(a)  such  that  Af{c  a(n))  is  PC  (i )  for 
every  n  >n0(a).  If  (A f  (c a(n ))  |  n  >p }  is  eventually  PC(i)  for  every  0<a<l, 
then  the  test  statistic  /  is  eventually  PC  (i).  □ 

It  is  now  natural  to  ask  whether,  for  a  monotone  test  statistic  / ,  the  properties 
“f  is  eventually  PC  (i)”  and  “/  is  SSC  (/)”  are  equivalent  Again  this  is  not  true 
in  general,  as  shown  by  the  behavior  of  f6  and  /7  above,  but  it  is  more  nearly  true: 
although  fs  is  not  PC(1)  for  some  values  of  n,  a,  it  is  both  eventually  PC(1)  and 
SSC(l).  It  is  interesting  to  note  that  this  is  indeed  true  for  all  admissible  test 
statistics  f .  More  precisely,  we  introduce  the  following  definition: 

DEFINITION  4.2.  For  fixed  a,  the  sequence  {A^(ca(n))  |  n  >p}  is  said  to  be 
eventually  admissible  for  the  original  testing  problem  (1.5)  if  there  exists  nfa) 
such  that  Ay(c  a(n ))  is  an  admissible  acceptance  region  for  (1.5)  whenever 
n  >n  !(a).  If  (A y(c  a(n ))  |  n  >p  }  is  eventually  admissible  for  every  0  <  a  <  1,  then 
the  test  statistic  /  is  eventually  admissible  for  the  testing  problem  (1.5).  □ 
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The  proof  of  the  following  theorem  is  based  on  Schwartz’s  (1967b)  necessary 
condition  for  the  admissibility  of  an  invariant  acceptance  region  for  problem 
(1.5) — see  Anderson  and  Perlman  (1988). 

THEOREM  4.3.  Let  /  be  a  monotone  invariant  test  statistic.  If  /  is  eventually 
admissible  for  the  testing  problem  (1.5),  then/  is  SSC(l)  and  eventually  PC(1).  □ 

Schwartz’s  (1967b)  sufficient  condition  for  admissibility  implies  that  the  level 
a  tests  based  on  /i,/2>/3,/s  are  admissible  for  every  n  and  a,  hence  a  fortiori 
fi’fi’fi’fs  are  eventually  admissible.  It  has  already  been  noted  that 
fi<  fi’f'i >  fs  are  SSC(l)  and  eventually  PC(1),  in  agreement  with  Theorem  4.3. 
Both  /4  and  /7  are  eventually  inadmissible  (although  admissible  for  sufficiently 
small  n  or  a)  and  both  are  PC(1),  but/4  is  SSC(l)  whereas  /7  is  not  SSC(l).  The 
level  a  tests  based  on /6,/g,  and/9  are  inadmissible  for  every  n  and  a  and  none  of 
y*6>  /s>  f 9  is  eventually  PC(1),  but  f6  is  SSC(l)  whereas  /8,/9  arenotSSC(l). 

Thus,  neither  the  requirement  that/  be  eventually  PC(1)  nor  the  requirement 
that  /  be  SSC(l)  distinguishes  among  invariant  tests.  For  fixed  n  and  a,  however, 
the  requirement  that  /  be  PC(1)  for  n  and  a  is  nor  satisfied  by  every  admissible  test 
statistic.  In  particular,  the  Bartlett-Nanda-Pillai  statistic  f5  fails  to  be  parameter 
consistent  unless  its  critical  value  satisfies  ca(p,r,n;f5)  <  1,  i.e.,  unless  n  or  a  is 
sufficiently  large.  (Values  of  c a(p  ,r  ,n;/s)  are  tabulated  in  Anderson  and  Perlman 
(1988).)  Despite  the  facts  that  the  Bartlett-Nanda-Pillai  test  is  admissible  (Schwartz 
(1967b)),  proper  Bayes  (Kiefer  and  Schwartz  (1965)),  locally  most  powerful 
invariant  and  locally  minimax  (Schwartz  (1967a)),  and  robust  (Olson  (1974)),  its 
failure  to  be  parameter  consistent  is  a  potentially  serious  drawback.  Injudicious  or 
routine  use  of  such  a  test  (for  example,  in  a  statistical  computer  package)  could 
result  in  failure  to  detect  a  sizable  departure  from  the  null  hypothesis  H0. 
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